(12) INTERNA! JONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(19) World Intellectual Property 
Organization 
Iniemalional Bureau 

(43) International Publication Date 
27 January 2005 (27.01.2005) 




PCT 



(10) International Publication Number 

WO 2005/008583 A3 



(51) International Patent Classilication^: G06T 7/20, 1 1/(M) 

(21) International Applicaticm Nujnber: 

PCT/IL20()4/(K)()632 

(22) International Filing Dale: 14 July 2004 (I4.()7.20()4) 

(25) Filing Language: English 

(26) Publication Language: English 



(30) Priority Data: 
60/488,382 



2 1 July 2003 (2 1 .07.2003) US 



(71) Applicant (for all designated States except US): PAIEON 
INC. rUS/USJ; 747 3rd Avenue, New York, New York 
10017-2803 (US). 

(72) Inventors; and 

(75) Inventors/Applicants (for US only): ZAKKH, Michael 



fn^yiL]; 9 Oranim Su^el, 54052 Givai Shmuel (IL). 
KLAIMAN, Moshe fn^ILI; 21/2 Marava Sirect, 70700 
Cjcdera (IL). 

(74) Agent: REINHOLD CORN AND PARTNERS; RO. 
Box 4060, 61040 Tel Aviv (IL). 

(81) Designated States (unless nthervi'ise indicated, for every 
kind of national protection available): AE, AG, AL, AM, 
AT, AU, AZ, BA, BB, BG, BR, BW, BY, BZ, CA, CH, CN, 
CO, CR, CU, CZ, DE, DK, DM, DZ, EC, EE, EG, ES, R, 
GB, GD, GE, GM, GM, TO, HU, ID, IL, IN, IS, JP, KE, 
KG, KP, KR, KZ, LC, LK, LR, LS, LT, LU, LV, MA, MD, 
MG, MK, MN, MW, MX, MZ, NA, NT, NO, NZ, OM, PG, 
PH, PL, PT, RO, RU, SC, SD, SE, SG, SK, SL, SY, TJ, TM, 
TN, TR, T\\ TZ, UA, UG, US, UZ, VC, VN, YU, ZA, ZM. 

zw. 

(84) Designated States (unless otherwise indicated, for every 
kind of regional protection available): ARIPO (BW, GH, 
GM, KE, LS, MW, MZ. NA, SD, SL, SZ, TZ, UG, ZM, 

[Continued on next page] 



(54) Title: METHOD AND SYSTEM FOR IDENTIFYING AN OPllMAL IMAGE WITl-nN A SERIES OF IMAGES TOAT 
DEPICn* A MOVING ORGAN 



INPUT: Series of images depicting a moving organ. Motion is | 
comprised of a cyclic motion and olher, non-cyclic, motion I 



Compute overall motion forwll the sequence of images 




"1 



Compute cycle period from 
ovcraH motion duta (for e.\amp!e 
by specirai analysis) 



Compute non-cyclic moiion 



(57) Abstract: A method and system for quantifying a cyclic mo- 
tion within a series of images depicting a moving object subject to 
composite moiion containing a cyclic component and a non-cyclic 
component of lower frequency than the cyclic component. Com- 
posite moiion is computed. The non-cyclic component is com- 
puted as the integral of moiion over a moiion cycle. The non- 
cyclic component is subtracted from Ihe composite moiion so as 
to obtain the cyclic component. 



BEST MAILABLE COPY 



< 

00 
IT) 



in 
O 



Subiraci nun-cyclic motion from overall 
moiion to obtain cyclic motion 






Use mohon values 
(cyclic, non-cvclic 
nnd overall) for direct 
measuremenrs 


1 


i 


Define an cvenlNs ihai is motion function 
related (miu, max, ...) 






Find candidates (per event) 
that match eveni\s 



Find the instance (per 
event) thai matches the 
even i\s 



T 



Find one instance (per 
event) from candidates as 
per heuristic decision 



The imjige (per event) thai correlaies lo this instance 
is\are the optimal imageXs for processing 



_2005008583A3„L> 



INTERNATIONAL SEARCH REPORT 



tnternanvnai Application No 

PCT/IL2004/000532 



A. CLASSIFICATION OF SUBJECT MATTER 

IPC 7 G06T7/20 



G06T11/00 



According lo Iniemaiional Paten! Classlficailon (IPC) or lo boih national classiricallon and IPC 



B. RELDS SEARCHED 



Minimum documentation searched (classificalion system lollowed by classification symbols) 

IPC 7 G06T GOIR 



Documeniation searched other than minimum documenlation to the exienl thai swch documents am included in the fields searched 



Elecironic data base con5aj»ed during the iniemaiional search (name ol data base and. where practical, search terms used) 

EPO-Internal. WPI Data, PAJ. INSPEC, COMPENDEX, IBM-TDB, BIOSIS, EMBASE 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category 



Cbalion of document, with indication, where appropriate, o( Ihe relevant passages 



Relevant to claim No. 



EP 0 885 594 A (MEDINOL LTD) 
23 December 1998 (1998-12-23) 

abstract 

page 2, line 42 - line 50 

page 4, line 20 - page 5, line 23 

page 7, line 30 - line 45 

page 10, line 38 - line 44 

page 13. line 22 - page 15, line 14 

figures 2,9 

-/— 



1,2. 

5-10,12, 
13.16-23 

3,4.11. 
14,15 



X I Further documents are listed in Ihe continuation of box C. 



iO 



Patent family members are listed In annex. 



Special categories of cited documents : 

'A' document defining the general stale of Ihe art which is not 
considered lo be of particular relevance 

E' earlier document but pubtlstied on or afterthe international 
filing date 

'L* document which may throw doubts on priority daim(s) or 
which is ciied to establish me publication date of another 
citation or other special reason (as specified) 

documant referring to an oral disdosure, use. exhibition or 
other means 

*P* document published prior to the iniemaiional filina date but 
later than the priority date claimed 



•T* later document published after the international filing date 
or priority date and not in conflict with the application but 
cited lo understand Ihe principle or theory underlying the 
invention 

•X' document of particular relevance; the claimed invention 
cannot be considered novel or cannot be considered to 
involve an inventive step when Ihe docurr^ent ts taken alone 

*Y* docunnent ot particular relevance; the claimed invention 

cannot be considered to involve an inventive step when the 
document is combined with one or more other such docu- 
ments, such combination being obvious to a person skilled 
in the art. 

*&* document member of the same patent family 



Date of the actual completion of the intemational search 



30 September 2004 



Date of mailing of the intemational search report 



3 I 01. 2005 



Name and maiBng address of the ISA 

European Patent Oti'ice, P.B. 5818 Patent taan 2 
NL - 2280 HV Rl|swijk 
Tel. (♦31-70) 340-2040. Tx. 31 651 epo nt, 
Fax: (-f 31 -70) 340-3016 



Authorized officer 



Engels, A 



Form FCT/iSA/210 (sscono sheet) (January 2004) 



DOCID: <WO 2005008583A3_I_> 



International Application No. PCT/ I L20O4/ D00632 



FURTHER INFORMATION CONTINUED FROM PCT/I5A/ 21G 



This International Searching Authority found multiple (groups of) 
inventions in this international application, as follows: 

1. claims: 1-23 

processing a series of images depicting a moving object 
subject to composite motion for extracting a cyclic motion 
component; 



2. claims: 24-28 

analysis of a graphical representation of a cyclic motion 
for identifying those images in. a series of images that 
depict an event associated with cyclic motion of a moving 
obj ect ; 



DOCID; <WO 2005O08583A3J_> 



INTERNATIONAL SEARCH REPORT 

Information on patent family members 



Patent document 
cited in search report 



Publication 
date 



Internanvnal Application No 

PCT/IL2004/000632 



Patent family 
member(s) 



EP 0885594 



23-12-1998 



Publication 
date 



US 


5924990 A 


20-07-1999 


AT 


283497 T 


15-12-2004 


AT 


236574 T 


15-04-2003 


AU 


757314 82 


13-02-2003 


AU 


7307498 A 


14-01-1999 


AU 


8032298 A 


04-01-1999 


BR 


9814778 A 


26-06-2001 


CA 


2240651 Al 


19-12-1998 


CZ 


9801901 A3 


13-10-1999 


DE 


19827460 Al 


24-12-1998 


DE 


69813087 Dl 


15-05-2003 


DE 


69813087 T2 


05-02-2004 


DE 


69827857 Dl 


30-12-2004 


DK 


885594 T3 


14-07-2003 


EE 


9800196 A 


15-02-1999 


EP 


1227342 Al 


31-07-2002 


EP 


0885594 A2 


23-12-1998 


ES 


2192290 T3 


01-10-2003 


GB 


2326479 A ,B 


23-12-1998 


WO 


9857580 Al 


23-12-1998 


XL 


125000 A 


12-05-2004 


JP 


11151246 A 


08-06-1999 


NO 


982817 A 


21-12-1998 


PL 


326831 Al 


21-12-1998 


FT 


885594 T 


29-08-2003 


RU 


2238041 C2 


20-10-2004 


SG 


68672 Al 


16-11-1999 


SK 


87098 A3 


14-02-2000 


US 


6152878 A 


28-11-2000 


US 


6095976 A 


01-08-2000 



Form PCT/ISA/ZIO (patert tamily annex) (January 2004) 



KDOCID: <WO 2005008583A3_L> 



(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 

lilllllillillllllilillllllUllllllllllilllili^ 



(19) World Intelleciual Property 
Organization 
Iniemaiional Bureau 



(43) Internationa] Publication Date 
27 January 2005 (27.01.2005) 




PCT 



(10) International Publicalion Number 

wo 2005/008583 A2 



(51) International Patent Classilicalion^: 
(21) international Application Number: 



G06T 1/00 



PCT/IL2004/00()r)32 

(22) International Filing Dale: 14 July 2004 (14.07.2004) 

(25) Filing Language: English 

(26) Publication Language: English 



(30) Priority Data: 

60/488,382 



2 1 July 2003 (2 1 .07.2003) US 



(71) Applicant {for all designaied Stairs except US): PAIEON 
INC. f US/US J; 747 3rd Avenue, New York, New York 
10017-2803 (US). 

(72) Inventors; and 

(75) Inventors/Applicants (for US only): ZARKH, Michael 



IFL/IL]; 9 Oranim Street, 54052 Givai Shmuel (IL). 
KLAIMAN, Moshe flL/IL]; 21/2 Ilarava Street, 70700 
Gcdcra (IL). 

(74) Agent: REFNHOLD COHN AND PARTNERS; P.O. 

Box 4060, 61040 Tel Aviv (IL). 

(81) Designated States (unless oiheryi'ise indicated, for every 
kind of national protection available): AE, AG, AL, AM, 
AT, AU, AZ, BA, BB, BG, BR, BW, BY, BZ, CA, CII, CN, 
CO, CR, CU, CZ, DE, DK, DM, DZ, EC, EE, EG, ES. M, 
GB, GD, GE, GH, GM, Hl^, IRJ, ID, IL, IN, IS, JP, KE, 
KG, KP, KR, KZ, LC, LK, LR, LS, LT LU, LV, MA, MD, 
MG, MK, MN, MW, MX, MZ, NA, Nl, NO, NZ, OM, PG, 
PH, PL, IT, RO, RU, SC, SO, SE, SG, SK, SL, SY TJ, TM, 
TN, TR, TT, TZ, UA, UG, US, UZ, VC, VN, YLl, ZA, ZM, 
ZW. 

(84) Designated States (unless otherwise indicated, for every 
kind of regional protection available): ARIPO (BW, GH, 
GM, KE, LS, MW, MZ, NA, SO, SL, SZ, TZ, UG, ZM, 

(Continued on next page} 



(54) Title: METHOD AND SYSTEM FOR IDENTIFYING OPTIMAL IMAGE WITHIN A SERIES OF IMAGES THAT DEPICT 
A MOVING ORGAN 



IKPUT: Scries of tniayes depicting a moving organ. Motion is 
comprised of a cyclic motion and other, non-cyclic, motion 




Yes 



Compute cycle period from 
overall motion data (for example 
by spectral analysis) 



Compute non-cyclic motion 



(57) Abstract: A method and system for quantifying a cyclic mo- 
tion within a series of images depicting a moving object subject lo 
composite motion containing a cyclic component and a non-cyclic 
component of lower frequency than the cyclic component. Com- 
posite motion is computed as well as the non-cyclic component as 
the iniegral of motion over a motion cycle. The non-cyclic com- 
ponent is subtracted from the composite motion so as to obtain the 
cyclic component. 



< 



r4 

o 



Subtract nou-cyclic motion from overall 
motion to obtain cyclic motion 








Use motion values 
(cyclic, non-cyclic 
and overall) for direct 
measurements 


I ' 


1 




Define an cvcnt'\s that is motion function 
related (min. max, ...) 




Find candiclyres (per tvcutl 
tbat inatcli eveiu'\s 



Find the instance (per 
event) that matches the 
evem\s 



Find one instance 0:>er 
event) from candidates as 
per heuristic decision 



The image (per event) ibai correlates to ilii5 instance 
I is\arc the optima) imageVs for processing 



DOCtD: <WO 2005008583A2_L> 



wo 2005/(K)8583 



PCT/IL20()4/(>00ri32 



-1- 

Method and system for identifying optimal image within a series of images that 

depict a moving organ 

FIELD OF THE INVENTION 

This invention relates to medical image processing devices. 

BACKGROUND OF THE INVENTION 

Medical imaging devices are often used to image moving organs. Cardiac image 
5 processing devices, in particular, are always used to image moving organs, either the heart 
(^da ultrasound imaging for example), or the coronaries (via angiography for example). 
Many of these imaging processing devices are used to quantify the motion either as an 
indication by itself or as pait of an image-processing algorithm. 

An image processing device for Left Ventricle Analj^sis is used to evaluate 

10 Ejection Fraction, which is the percentage of the blood pumped out during each heartbeat. 
Left Ventricle Analysis involves computing the Left Ventricle volume from an angiogram 
(taken from a cme-angio sequence of images). The Left Ventricle volume is computed 
once for die heait in its systoUc phase and once for the heart in its diastolic phase. Ejection 
Fraction is estimated from the ratio of these volumes. Identifying the systoHc and diastoUc 

1 5 images is part of the LVA procedure. 

Myocardium thicloiess and Heart Wall Motion are evaluated from Ultrasound 
Images to indicate heart failure conditions. Both procedures, again, involve the 
identification of systolic and diastolic instances. Furthermore, quantifying the object's 
motion could be directly used for Wall Motion evaluation. 

20 hitra- Vascular Ulti'asound (TVLTS) is a method of evaluating and anafyzing 

coronary defects by means of inserting an intra-vascular ultrasound device and imaging the 
vessel. rVUS measui*ements include measurements of the luminal ^^essel area. Estimation 
of the luminal area verj' much depends on the heart phase and results vary for different 
images depicting dijfferent stages within the cardiac cycle. Again, it is usefiil to identify the 

25 diastolic - or the nunimal movement instance - in order to perfomi measurements on the 
optimal image. 
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SUMMARY OF THE EWENTION 

It is an object of the invention to provide a method and system for imaging 
moving organs and to quantify the organ movements. 

The present invention provides several methods and sj'Stems that relate to the 
5 evaluation of an orgjan motion. 

According to tlie invention, there is provided a method for quantifjdng a cyclic 
motion within a series of images depicting a moving object subject to composite motion 
containing a cyclic component and a non-c5'clic component of lower frequency than the 
cyclic component, the method comprising: 
10 (a) computing the composite motion; 

(b) computing the non-cyclic component as the integral of motion over a motion 
cycle; and 

(c) subtracting the non-cyclic component from the composite motion so as to 
obtain the cyclic component 

15 The invention provides a novel method of evaluating the cyclic motion of an 

organ from a series of images of any source and provides implementations of such a 
method that decrease or eliminate motion artifacts. Specifically, we present a novel 
method and system for selecting optimal images for the process of 3D reconstruction of 
the coronaries. We fiirther provide a method and system for replacing the need for ECG 

20 Gating b)'^ an analysis of the heart movement 

A metliod for estimating the motion of an organ of a series of images comprises 
the folloAving operations. A medical imaging device acquires a series of images presenting 
an organ that is in motion. The motion is either of the organ changing shape (eventually, in 
a cyclic manner, regaining its original shape) or additionally of the organ changing 

25 location mthin the image (in angiography, for example, it is very common to move the 
patient's bed while imaging; as a result shifting the coronaries' location within tlie image). 
If a non-cyclic component is superimposed on the cyclic motion, the series of images are 
analyzed to separate the cyclic motion from the non-cyclic motion. Once these two types 
of motion are separated, the cyclic motion can be quantified. The quantified motion can 

30 now be used for dkect measurements or can be investigated to identify different events 
within the motion cycle. In some implementations, this investigation will point to an 
unage that is optimal in the sense that it represents minimal motion and thus j'ields 
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DETABLED DESCRIPTION OF EXEMPLARY EMBODIMENTS 
General principles 

The invention will be described with particular reference to the determination of 
an optunal image within a series of images depicting heart motion possibly contauiing 
5 "noise" caused, for example, by shifting of the operating table on which a patient is 
disposed. Before doing so, some general algorithms will JBrst be described. 

Fig. 2 is a flow chart showing the principal operations carried out by a method 
according to the invention for identifying an optimal image from a series of images that 
depict a moving object. 

10 Fig. 3 is a flow chart shov^dng a preferred embodiment of such a method for 

identifying optimal image from a series of coronary cine-angiography images; this image 
being the input for a tliree-dimensional reconstruction of a the coronaries. 

A series of images is received as input to from any imaging source. These images 
depict a mo\dng organ subjected to two types of motions. The first is a cyclic motion of 

15 the object itself, meaning that, within a certain time frame, the object restores its original 
shape and position. The second is the motion of the object within the scene (image), 
meaning that the object changes position due to change in the imaging position. In the 
preferred embodiment of this method shown in Fig. 3, the images are a series of coronary 
angiographic images, obtained during a catheterization procedure. The cyclic motion is the 

20 heartbeat and the second motion could be, for example, movement of the C-ARM table, 
causing a shift of the imaged coronary vessels in the image. 

First, the overall motion is computed for all sequence of images, using any 
method, for example optical flov/ or phase correlation. 

If the cyclic period is unlaiown, then the cycle period is computed from the 

25 overall motion data. One method of doing so is by spectral anafysis. The non-cyclic 
motion is computed, using overall motion and known or computed cycle period. A 
preferred embodiment of a known cycle period is the period of a cardiac cycle extracted 
from analysis of the ECG signal. 

The non-cyclic motion is subtracted from the overall motion to obtain the cyclic 

30 motion. In the preferred embodiment shown in Fig. 3, the heartbeat motion is obtained by 
subtracting the non-cyclic motion (mainly attributed to movement of the patient's bed) 
from the overall motion. 

.2005008583 A2_L> 
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this sequence are attributed to the organ's cyclic motion, but ai'e also attributed to otlier 
factors. 

If only cyclic motion is present, the first and the last images in this sequence must 
be identical, IMi=lMni+i. 
5 Difference between images representing composite motion can be computed, as 

known in the art, by, for example, optical flow or by applying phase correlation 
computation to pairs of successive images IMj and IMj+i, ie{l..m}. The result of this 
computation (for example the result of phase con-elation) is described dXj ,dYi and pi, 
where dXj and dYj are the shift between images (assuming a substantial part of the same 
10 pattern is present in both images) in X and Y axes respectively and pi is the correlation 
gi*ade. Pi may be used to enhance the flirtlier described algorithms. 
Let us define and compute the motion integration as: 









+ 












dY; 



15 meaning that in tlte first image, the motion integral is equal to zero. The motion integral 

for image i+l is equal to the motion integral for image / plus shift between images / and 

7+7, as computed b}^ tlie phase correlation. 

It is mathematically understandable that the integration of a cycUc motion, from 

image IMi to IMm+u is zero - if an object starts and ends in the same position, then the 
20 integration of the object movements (on X and Y axis) is zero (as shoAra in Fig. 1). Thus, 

if only oycYic motion is present, (Xm-M,Ym+i)=(Xi,Yi)=(0,0). 

Let (Xnc. Ync) be the integral of the non-cyclic motion, 

(XnCjYnc) = (Xm+uYm+l), 

25 This means that, given that tlie integral of cyclic motion is 0, (Xm+uYm+i) 

represent the residual motion that is attributed to non-cyclic movement. 

Assuming non-cyclic motion is consistent, or at the least that its frequency is 
lower than the cyclic motion frequency, we can subtract this motion from the overair 
30 motion: 

(Xi*,Yi*)=(Xi,Yi)-(XNc,YNc)*(i-l )/m, i=l ,2,. . .,m+l . 
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Dlm,e = max {Di,E} for all /=L.m+L 

We can relieve tlie requirement for knowing in advance an approximation for the 
least motion image by means of direct computation. If indeed the least motion instance is 
5 not unique, we can use heuristic criteria for distinction. For example, within a sequence of 
angiograms, depicting a cardiac cycle, it is easy to distinguish between the end-systole 
instance and tlie end-diastole instance, both representuig least motion, since the end- 
diastole instance is identified by presentmg the coronaiies in maximal spreading, while the 
end-systole instance is identified by presenting the coronaries in minimal spreading. 

10 

Preferred Embodiment. 

We suggest a preferred embodunent for an application of three-dimensional 
reconstruction of coronary vessels from a procedure of conventional angiography. In order 
to reconstruct a thi*ee-dimensional image of the arteries, it is necessary to obtain at least 

15 two two-dimensional images of the arteries m tlie same phase of the heartbeat, for 
example at end-diastole. Therefore, image acquisition is usually synchronized to an ECG 
signal. This procedure involves simultaneous recordings of tlie video signal firom the X- 
ray camera and the patient's ECG signal. We present here a novel method for identifying 
the end-diastole instance, equivalent to ECG-gating, without relying solely, if at all, on the 

20 ECG signal. 

Let IMj, IM2,.. IMn be n images of a catheterization-acquired run. 

Let m be the number of frames per cardiac cycle, either known in advance or 
computed as detailed in the above-described metliod for estimating the organ's motion. 
25 Let IMk be the approximate location of end-diastolic frame within the cycle, 

either IcnoAvn in advance or heuristically identified as detailed in the above-described 
method for obtaining the least motion image. 

EN'Ik-m/i^ IMk.in/2+1, IM k-ni/2+1 • - - IMk.m/2+m fotm a fiill cardiac cycle (for the sake of 
30 sunpUcity, let us presume that m is an even number). Differences between frames in this 
sequence are attributed to heart motion, but are also attributed to bed motion, iodine 
propagation and several other reasons. If only heait motion were present, the fii'st and the 
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Let: 

We can determine the end-systole point S, wliich is the one most distant from the 
approximated end-diastole point, meaning: 
5 Ds,m/2+i = max {Di.m/2+1} 

Minimum motion point - end-diastole, ED • is determined as most distant from 
systole point: 

Ds,ED = max {Ds,j}: 

1^ Selecting the IMed image per sequence of cine-angio images for the process of 

thi-ee-dimensional reconstruction will provide the optimal result, in terms of accuracy and 
precision, for the reconstruction and for vessel analysis. 

In the method claims that follow, alphabetic characters and Roman numerals 
used to designate claim steps are provided for convenience only and do not imply any 

1 5 particular order of perfomiing the steps. 

It wall also be xmderstood that the system according to the invention may be a 
suitably programmed computer. Likewise, the invention contemplates a computer progi'am 
being readable by a computer for executing the method of the invention. Tlie invention 
ftirther contemplates a machine-readable memoiy tangibly embodying a program of 

20 instructions executable by the machine for executing the metliod of the invention. 
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11. The metliod according to claim 10, wherein the selected image is closest to a 
predetermined approximation. 

12. The method according to claim 10 or 11, wherem the event is least motion. 

13. The method according to claim 12, for selecting angiograpliic images to 
5 participate in thi'ee-dimensional reconstiiiction of coronary vessels. 

14. The method according to claim 13, including deriving cycle period and approxi- 
mation for least-motion image from an analysis of an ECG signal. 

15. The method according to claim 13 or 14, including distinguisliing the end- 
diastole instance from the end-systole instance by the state of coronary vessel - maximal 

i 0 spreading versus minimal spreading, respectively. 

16. The method according to any one of claims 5 to 15 when used for selecting 
opthnal image or images for QCA analj^sis. 

17. The metliod according to any one of claims 5 to 15 when used for selecting 
optimal image or images for IVUS analysis. 

15 18. The method according to any one of claims 5 to 15 when used for selecting 
optimal image or images for LVA analysis. 

19. The method according to any one of claims 5 to 15 when used for selecting 
optimal image or images for Wall Motion analysis. 

20. The metliod according to any one of claims 5 to 15 when used for CT 
20 reconstruction. 

21. The method according to any one of clauns 5 to 15 when used , for MRI 
reconstruction. 

22. The method according to any one of clauns 5 to 15 when used for PET 
reconstruction, 

25 23. A system for quantifying a cyclic motion within a series of images depicting a 
moving object subject to composite motion containing a cycUc component and a non- 
cj^clic component of lower fi-equency than the cyclic component, the system comprismg: 
a composite motion unit computing the composite motion. 
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INPUT: Series of angiograms depicting the coronary arteries. 
Motion is comprised of the heartbeat (cyclic) and other 
movements, especially movement of the 
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Compute non-cyclic motion 
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Subtract non-cyclic motion from overall 
motion to obtain cj^clic, heartbeat motion 
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We define the event of least arteries motion 
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least-motion (end-diastole 
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Identify'' end-diastole as the 

image depicting a most 
relaxed state of the arteries 



FIG. 3 



The image that depicts the instance of least-motion 
most-relaxed state of the arteries is the cpptimal iinage f©F 
three-dimensional reconstruction 
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